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We analyze the interaction-induced renormalization of single-particle excitations 
in the two-dimensional Hubbard model at weak coupling using the Wick-ordered 



D ' version of the functional renormalization group. The self energy is computed for real 

H 

"^ ■ frequencies by integrating a flow equation with renormalized two-particle interac- 



ts ■ tions. In the vicinity of hot spots, that is points where the Fermi surface intersects 

the umklapp surface, self energy effects beyond the usual quasi-particle renormaliza- 
tions and damping occur near instabilities of the normal, metallic phase. Strongly 
enhanced renormalized interactions between particles at different hot spots generate 

1^ \ a pronounced low-energy peak in the imaginary part of the self energy, leading to 

^ I a pseudogap-like double-peak structure in the spectral function for single-particle 

T-H ! 

^O ' excitations. 

O 



PACS: 71.10.Fd , 71.27.+a , 71.10.Hf 



^ ■ I. INTRODUCTION 

o' 
o 



The most popular prototype model for the dynamics of electrons in the copper oxide 
planes of cuprate superconductors is the two-dimensional Hubbard model,— which captures 
two of the most prominent and universal features of these compounds: to be an antiferro- 
magnetic Mott insulator at half filling and a d-wave superconductor at sufficiently strong 
doping. It is very difficult to compute the low-energy properties of the two-dimensional 
Hubbard model at strong coupling, as relevant to the cuprates. In fact, close to half filling 
the model exhibits complicated non-trivial behavior already at weak coupling, due to a 
mutual interplay of singularities in the particle-particle and particle-hole channel, which 



cannot be captured by simple resummations of Feynman diagrams. 

A systematic analysis of this interplay has been initiated in recent years with the 
development of functional renormalization group (fRG) methods for interacting Fermi 
systems.— The starting point of this approach is an exact hierarchy of differential flow 
equations for Green or vertex functions, which is obtained by differentiating the generating 
functional of the system with respect to an infrared cutoff A or some other energy scale. 
Approximations are then constructed by truncating this hierarchy and parametrizing the 
vertex functions with a managable set of variables or functions. Already from the flrst fRG 
studies of the two-dimensional Hubbard model, which focussed on the flow of the two- 
particle vertex F and various susceptibilities in the one-loop approximation, a number 
of interesting conclusions could be drawn.i^^ In particular, the expected existence of 
superconductivity with d^^-y^ symmetry was flrmly established, at least at weak coupling. 
More importantly, it became clear that close to half fllling correlations in the particle- 
particle and particle-hole channel are both strong and mutually influence each other. 

The effect of correlations on single-particle excitations is seen most clearly in the self 
energy S(k, cu), for frequencies u) on the real axis. Zanchi^ computed the flow of the 
frequency derivative of the self energy on the Fermi surface, which determines the wave 
function renormalization factor Z = [1 — 9^ReS(ki?,a;)|^^Q]~^, for the special case of 
a half-fllled Hubbard model with pure nearest neighbor hopping. In Fermi liquids Z 
yields the weight of the quasi-particle peak in the spectral function. Zanchi found that 
the Z-factor is strongly reduced when the renormalized vertex increases, the suppression 
being strongest near the van Hove points. Subsequently, Honerkamp and Salmhofer— 
computed the flow of the Z-factor for the two-dimensional Hubbard model at as well as 
away from half fllling. While conflrming an anisotropic reduction of Z, they observed 
that this suppression emerges much slower than the divergence of dominant interactions. 
In addition, they pointed out that the self energy does not signiflcantly alter the vertex 
renormalization in the perturbative regime, that is before the interactions grow strongly. 
An fRG calculation of the imaginary part of the self energy at zero frequency, ImS(ki7, 0), 
was performed by Honerkamp.— In Fermi liquids this quantity is directly related to the 
decay rate of quasi-particles. For a hole-doped Hubbard model with a concave Fermi 



surface that almost touches the van Hove points he found that IraUCkF, 0) is moderately 
anisotropic, with larger values for momenta closer the van Hove points. A closer look 
reveals that the imaginary part of the self energy is actually largest near a hot spot, 
foreshadowing already a main feature of our results presented below. The first fRG 
calculation of the full frequency dependence of the self energy has been performed very 
recently by Katanin and Kampf.-' They computed the real and imaginary parts of S(k^, uj) 
and the resulting spectral function ACkp, u) for the Hubbard model with nearest and next- 
to-nearest neighbor hopping at temperatures above T* and observed the destruction of 
Fermi liquid behavior near T* due to structures indicating the formation of a pseudogap, 
in particular for kp close to a van Hove point. For kp at short but finite distance from a 
van Hove point they found an additional peak between the two maxima that confine the 
pseudogap in the spectral function. Moving even further away from the van Hove point 
toward the diagonal of the Brillouin zone the three peaks gradually merge into one single 
peak.— 

In this work we present a calculation of the full frequency dependence of the self energy 
and the resulting spectral function for the two-dimensional Hubbard model as obtained 
from the Wick- ordered version of the fRG.-^ This scheme has the advantage that only low- 
energy states contribute to the flow at low-energy scales, which leads to a better control 
of the necessary Fermi surface projections of momentum dependences compared to other 
fRG versions.i^*^ We do the calculation for frequencies directly on the real axis, thus 
avoiding analytical continuation. In agreement with previous studies we obtain marked 
deviations from Fermi liquid behavior for Fermi momenta near van Hove points, and also 
near other hot spots in the case of filling factors above the van Hove singularity: the 
imaginary part of the self energy develops a pronounced peak at low frequencies, which 
leads to a double-peak structure in the spectral function reminiscent of a pseudogap. This 
double peak gradually transforms into a single peak when moving away from the hot spot 
or van Hove regions, and also upon raising the temperature. Unlike Ref.|9|we never obtain 
a third peak in the pseudogap valley. 



II. MODEL 

We consider the one-band Hubbard model defined by the Hamiltonian 

with a local repulsion f/ > and hopping amplitudes tjj/ = —t between nearest neighbors 
and tjj' = —t' between next-to-nearest neighbors on a square lattice. The corresponding 
dispersion relation 

Ck = — 2t(cos kx + cos ky) — At' cos k^ cos ky (2) 

has saddle points at k = (0,7r) and (7r,0), leading to logarithmic van Hove singularities 
in the non-interacting density of states at energy evH = 4t'. 

III. METHOD 

We use the Wick-ordered version of the fRG as proposed by Salmhofer,— which has been 
concisely described and used in the context of an explicit one-loop computation of the 
two-particle vertex for the two-dimensional Hubbard model.'^ In this scheme, to second 
order in the two-particle vertex the flow of the self energy E^ is determined by a single 
two- loop Feynman diagram, shown in Fig. 1, and the flow of the two-particle vertex F^ is 
obtained to second order by evaluating the one-loop diagram in Fig. 1. The internal lines 
without slash in the Feynman diagrams correspond to the bare propagator 

ikQ - ^k 
where ^k = Ck — /U and A > is the cutoff; the slashed hues represent d^D^, which is 
proportional to 5(A — |^k|)- Note that D^ has support only on a shell of energetic width 
A around the Fermi surface, that is for momenta k with |^k| below scale A. A first order 
contribution to S^ (Hartree term) has been absorbed in a shift of the chemical potential, 
such that the leading contributions are quadratic in F^. We note that in the full hierarchy 
of flow equations of the Wick-ordered fRG also one-particle reducible terms contribute. 



However, these have to be considered only if the flow is computed beyond second order 
inT'^. 

For a spin-rot at ionally invariant system the vertex is uniquely determined by two 
functions depending on momenta only, namely a singlet component F^ and a triplet 



component F^, as described in Ref. 



n' 

• i 



In the present work the vertex is defined with 

The flow equation for the 



an extra factor 1/2 compared to the conventions in Ref. 
vertex can only be solved if the momentum and frequency dependence of the vertex 
is parametrized by a reduced set of variables. The frequency dependence of F^ can be 
neglected without much damage because it is absent in the bare interaction, and irrelevant 
in the sense of power counting in the limit A — > 0. The frequency integrals in the Feynman 
diagrams in Fig. 1 can then be carried out analytically. The resulting one-loop flow 
equations for F^(k']^, k2; ki, k2) and Ff (k'^, k2; ki, k2) given in Ref. i^ are supplemented by 
the flow equation for the imaginary part of S^(k, u) for arbitrary real frequencies u given 
as 

X {[F^(k,q;p,k + q-p)]2 + 3[Ff(k,q;p,k + q-p)]2} 

X F(q;p,k+q-p)5(c.;-^p + ^q-^k+q_p) (4) 

with the combination of Bose (6) and Fermi (/) functions 

^(q, p, p ) = [fiQ - /(^p')] [/(^p) + K^q - ^p')] • (5) 

In contrast to its appearance F is symmetric in p and p'. The real part of S^ can be 
determined from ImS^ via a Kramers- Kronig relation. The momentum dependence of the 
vertex is discretized similarly to previous fRG calculations^^ by dividing the Brillouin 
zone into patches as shown in Fig. 2. The shape of the patches accounts for the fact 
that the momentum dependence normal to the Fermi surface is irrelevant in the sense of 
power counting, in contrast to the tangential momentum dependence. The flow equations 
are integrated from A = Aq = maxk |^k| to A = 0. The initial conditions are F^" = U, 
Fj ° = 0, and S^° = 0. From the result for the self energy S^ for A ^ we compute the 
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spectral function A{h,uj) = — 2ImG*(k, a;). To avoid the computationally expensive self- 
consistent determination of the interacting Fermi surface, we remove Fermi surface shifts 
by subtracting ReS(k^, 0) from the self energy in the Dyson equation G^^ = Gq^ — S . 

IV. RESULTS 

We focus on the two-dimensional Hubbard model with a negative next-to-nearest neigh- 
bor hopping amplitude t' = — O.lt and two different average densities, n = 0.92 and n = 1, 
respectively. In both cases the Fermi surface given by ,^k = is closed around (tt, tt) and 
has finite curvature everywhere, see Fig. 2. In the following the so-called hot spots, that 
is the points where the Fermi surface intersects the umklapp surface, will play a special 
role. The momenta on the umklapp surface obey the relation ek+q = Ck with Q = (tt, tt). 
For n = 0.92 the Fermi surface almost touches the van Hove points, and the hot spots are 
very close to the latter. By contrast, for n = 1 the hot spots are well separated from the 
van Hove points. We set t = 1, that is all energies are given in units of t. We fix a low 
finite temperature T = 0.05 and choose the bare interaction U such that the couplings 
fiow to strong values but do not diverge for A ^ 0. In other words, the temperature T is 
slightly above the scale T* at which the one-loop fiow diverges. 

The complete frequency dependence of the imaginary part of the self energy ImS(ki?, uj) 
for high and low energies is shown in Fig. 3, for n = 0.92 and a Fermi momentum 
k^ which coincides with a hot spot. The fRG result is compared to the result from 
standard second order perturbation theory, which is obtained from the fiow equation (4) 
when the renormalized vertex F^ is replaced by the bare interaction. While the vertex 
renormalization has virtually no infiuence on the self energy at high frequencies, it leads 
to a pronounced peak in ImS at low frequencies, which is absent in perturbation theory. 

We now focus on the low-energy behavior. In Fig. 4 we plot ^(k^,^;) at small fre- 
quencies for n = 0.92 and six different Fermi momenta k^ situated in patches 1, 2, . . . , 6, 
respectively. We show two different fRG results: in one case the bare interaction U has 
been tuned such that the renormalized vertex reaches a peak value F^'^^ = 15, in the other 
^max ~ ISO- ■'■f ^^^ peak in ImS(k^,cj) observed already in Fig. 3 is sufficiently strong. 



it causes a sign change in the slope of ReS(kir,u;) at small uj. For t}^ReS(ki?,u;)|^^Q > 1 
this leads to a double-peak structure reminiscent of a pseudogap in the spectral func- 
tion A(k^,a;), shown in Fig. 5. For a less strongly renormalized vertex, as well as for 
momenta away from the hot spot, the peak in ImS(k^,C(j) disappears gradually. As a 
consequence, the double peak in the spectral function is well developed only for ki? near 
the hot spot, and for a comparatively large renormalized vertex, that is for T rather close 
to T*. Otherwise, the self energy merely leads to an angle-dependent lifetime broadening 
of the quasi-particle peak in A{\lp^lS). 

For electron density n = 1 the hot spot is situated in the third patch (see Fig. 2). 
Results for S(kir, lS) and A(ki?, uS) are shown in Figs. 6 and 7, respectively. The fRG again 
yields a peak in ImS(k^,C(j) around c<j = for a hot spot momentum kp, but only for 
^max — 150. Moreover, its absolute value is much smaller compared to the case n = 0.92, 
and it is less pronounced. Consequently, the corresponding spectral function exhibits only 
a mild dip. For weaker interactions or kp away from the hot spot A(ki?, u) has a more or 
less conventional shape, with a single lifetime-broadened peak. This broadening is larger 
for momenta next to the hot spot, since the low-frequency part of ImS(ki?,a;) is still 
strongly enhanced, albeit not peaked, in these regions. 

Note that there is nothing special about the density n = 1 (half-filling) at weak coupling 
if t' 7^ 0, in contrast to the case of strong coupling, for which the half-filled Hubbard model 
becomes a Mott insulator, and in contrast to the special case with pure nearest neighbor 
hopping {f = 0), for which the Fermi surface coincides with the umklapp surface at half- 
filling. The fact that the anomaly we observe in the self energy is more pronounced for a 
Fermi surface close to van Hove filling reflects the relevance of the van Hove singularities 
in the weak coupling limit. 

To understand the mechanism behind the low-energy singularity in ImS(k, a;) for hot 
spot momenta khs, we consider the flow eq. (4) in the limit T ^ and uj ^ 0. At T = 
the function -F(q, p, p') reduces to a step function which constrains the relative signs of 



^q, ^p, and Cp'- 

-1 for ^q<0; ^p.^p' >0 

i"(q, P, P ) = < -1 for ^q > ; Cp, Cp' < (6) 

else 
v 

Energy conservation imposed by the (5-function in Eq. (4) then imphes that in the hmit 
C(j ^ the internal momenta q, p, p' = k + q — p are all constrained to the Fermi 
surface. If in addition k is a Fermi momentum, all ingoing and outgoing momenta of F^ 
in Eq. (4) have to lie on the Fermi surface. At low finite T or cj this constraint is slightly 
relaxed such that the internal momenta can be taken from a thin shell of energetic width 
of order T ox uj around the Fermi surface. Since one of the internal momenta p on the 
right hand side of the flow equation is forced to sit at scale |^p| = A by the A-derivative 
of the cutoff functions, the flow of ImS"^ receives sizable contributions only when A has 
reached the scale T ox uj. In two dimensions there are two generic types of momentum- 
conserving interaction processes with all ingoing (ki,k2) and outgoing (k']^,k2) momenta 
on the Fermi surface: forward scattering processes (k'^^ = ki or k'^ = k2) and Cooper 
processes (k2 = — ki). In lattice systems also certain umklapp processes in which a 
reciprocal lattice vector is added to the total momentum can fulflll this condition. They 
can play an important role, as discussed below. 

Looking at the renormalized interactions obtained from the one-loop flow of F"^ for our 
two choices of densities, n = 0.92 and n = 1, we flnd that in both cases the dominant 
effective interactions are in the spin singlet channel with ingoing and outgoing momenta 
near hot spots. Among these, interactions with a momentum transfer Q = (tt, tt) are 
leading. For n = 0.92 the vertex is maximal for the scattering of a Cooper pair with 
both partners on hot spots to two other hot spots with a momentum transfer Q, that is 
(khs, — khs) t-^ (k'j^s) "'^hs) with k'j^g — khs = Q. Combined Cooper and forward scattering 
processes (khs, —khs) ^-^ (khs, —khs) are also among the strongly enhanced interactions 
at that density. Computing the contributions to ImS from different scattering processes 
separately reveals that the peak at hot spots is indeed generated mostly by Cooper and 
Q-transferring interactions, where the latter are slightly more important. For n = 1 the 
leading interaction is an umklapp process, (khs, khs) ^-^ (khsi^hs) with kj^^ — khs = Q, 
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closely followed by a Q-transferring Cooper process on hot spots. However, the total 
contribution from Cooper processes (including the Q-transferring one) to the peak in 
ImS is very small compared to the total contribution from Q-transferring processes in that 
case. The importance of Q-transferring and Cooper interactions in the two-dimensional 
Hubbard model is well known from earlier fRG and other studies, and reflects strong 
antiferromagnetic and superconducting correlations, respectively. 

The above mechanism for the destruction of fermionic quasi-particles close to the Fermi 
surface is completely different from the mechanism leading to Luttinger liquid behavior 
in one-dimensional electron systems.— In one-dimensional (gapless) electron liquids the 
renormalized interactions remain finite, while the self-energy receives nevertheless en- 
hanced contributions at low energy due to a large phase space for scattering processes. 
However, linJ]{kF,uj) does not develop a peak at cj = 0, it just vanishes much slower for 
T — i> 0, C(j — > than in a Fermi liquid. At T = 0, one obtains lTaE{kF,uj) oc |l(j|^~'', where 
?7 is a positive exponent which is relatively small for systems with short-range interactions. 
The spectral function A{k,uj) exhibits a single peak centered at cj = for k = kp, and a 
split peak at finite energy for k away from the Fermi point.— This peak splitting is not 
related to pseudogap behavior, it is rather due to spin-charge separation. Applying our 
fRG scheme to a one-dimensional Luttinger liquid system at weak coupling yields finite 
renormalized interactions, due to a well-known cancellation of singularities in the particle- 
particle and particle-hole channel. For the self-energy as computed from the flow equation 
(4) in one dimension this implies a linear frequency dependence, ImS(A;p',co') oc \uj\, yield- 
ing a single broadened peak in A{kp, u) as expected, but missing the anomalous exponent 
rj. The latter could be captured by feeding the self-energy back on the right hand side 
of its flow equation, most simply via a wave function renormalization, as described for 
the one-particle irreducible fRG in Ref. lS- The situation is different if umklapp scattering 
(so-called (?3-terms) are relevant, as in the one-dimensional Hubbard model at half-filling. 
In that case renormalized interactions grow strong and a charge gap opens, leading to 
a gap in the spectral function A{kp,u!), which is clearly signalled by a suppression of 
spectral weight at low frequencies within the fRG scheme. 

Our results agree only partially with those obtained by Katanin and Kampf,— who 



10 

computed the full frequency dependence of the self energy from the one-particle irreducible 
version of the fRG,— with a truncation of the exact hierarchy of flow equations at the 
same order as in our calculation. For a Fermi momentum very close to a van Hove point 
they obtained a single peak in ImS resulting in a double peak in the spectral function, in 
agreement with our calculation. However, for Fermi momenta at some moderate distance 
from the van Hove point they found a double peak in ImS leading to a three-peak structure 
in the spectral function, which we do not observe in our results. In both approaches. 
Wick-ordered as well as one-particle irreducible fRG, the peaks in ImS are generated by 
strongly enhanced effective interactions, and are thus most pronounced at temperatures 
close to the instability scale T*. The different structure of ImS obtained from the two 
fRG versions can be traced back to a different division of momentum integrals in high 
and low energy modes. In both cases the energy variable uj is related to the excitation 
energies on internal lines: cu = ^p + ^p/ — ^q. In the Wick-ordered scheme all excitation 
energies are at or below the cutoff, that is |^| < A. As a consequence, the flow for 
ImS^(k, cij) with |ct;| > stops completely for A < |c(j|/3. In a situation where the vertex 
T^ grows strongly at very small scales, A ^ 0, the imaginary part of the self energy 
therefore receives the corresponding large contributions only for small frequencies, which 
leads necessarily to a single peak with maximal height at cu = 0. By contrast, in the 
two-loop diagram generating the flow of ImS^ in the one-particle irreducible scheme one 
excitation energy is at scale A, the second at a scale A' > A, and the third one above the 
scale A', where A' is the scale of the vertex F^ in the diagram. Note that the two- loop 
diagram for the self energy in the one-particle irreducible fRG is non-local in the cutoff. 
The crucial point is that not all excitation energies are bounded by the cutoff. Therefore, 
the flow of ImS^(k, uj) does not stop for A ^ |ci;|, and the strong enhancement of F^ for 
A' — > does not necessarily enhance ImS^(k, cj) most strongly for the smallest energies. 
However, to obtain a double peak in ImS with maxima at finite frequencies, say u^-, the 
contribution with at least one excitation energy of order uj± must be particularly large. 
This is possible only if the vertex is not enhanced most strongly for cases with all ingoing 
and outgoing momenta on the Fermi surface, but rather for some interaction processes 
with at least one momentum away from the Fermi surface. One should also keep in mind 
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that contributions from momenta away from the Fermi surface might be overestimated 
in the one-particle irreducible fRG with a momentum discretization by patches, since the 
dependence of the vertex T^ normal to the Fermi surface is weak only within a shell 
|.^k| < A in momentum space. 

We finally note that the Z-factor defined as Z = [1 — dujReT,(kF,uj)\^^Q]~^ does not 
provide sufficient information to parametrize the low-frequency dependence of the self 
energy, if the imaginary part develops a low-energy peak. Furthermore, Z does not 
provide any information on the low-energy weight in the spectral function in that case. 
Since the slope of ReS becomes positive near the hot spots, the Z-factor does not even 
satisfy the inequality < Z < 1. 

V. CONCLUSION 

In summary, we have presented a weak-coupling calculation of the self energy S(k, cj) 
for the two-dimensional Hubbard model with a renormalization of effective interactions 
which fully treats the interplay of particle-particle and particle-hole channels. Strong 
renormalizations of single-particle excitations via self energy effects were observed for 
momenta near hot spots on the Fermi surface, if the system is close to an instability 
where the renormalized interactions are strongly enhanced. In this case the quasi-particle 
peak in the spectral function A(k, uj) is replaced by a double-peak structure reminiscent 
of a pseudogap. Further more, the fRG flow enables us to identify the correlations which 
are responsible for this feature, and thus offers substantial insight into the underlying 
microscopic mechansim. Evidently, the perturbative truncation of the flow equations 
becomes more and more uncontrolled in the regime of large renormalized interactions. 
We can thus only conclude that single-particle excitations near hot spots are strongly 
affected by interactions, but higher order contributions to the self energy are expected to 
become important at a certain point. 

The tendency to pseudogap formation found in the above weak-coupling calculation 
is restricted to a narrow temperature scale and also to a narrow regime in momentum 
space, namely close to the hot spots. In contrast, the pseudogap phenomena in hole-doped 
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cuprate superconductorsi^ are observed in a much wider temperature range, and are seen 
almost everywhere on the Fermi surface. This is undoubtedly due to the fact that the 
interaction in these materials is strong. Indeed, numerical results for the two-dimensional 
Hubbard model at large U obtained from quantum cluster approaches^ provide evidence 
for extended pseudogap behavior near half-filling. 

Our observation that for moderate interaction strength the pseudogap formation is 
favored at hot spots is in accordance with recent results obtained by Kyung et al.— on 
the basis of the two-particle self-consistent (TPSC) approach, as well as with the find- 
ings by Senechal and Tremblayi^ from cluster perturbation theory. In both references 
the focus is put on electron-doped cuprate superconductors, and it is argued that these 
compounds may be more weakly coupled at optimal doping, as opposed to the hole-doped 
materials. Therefore, it seems more appropriate to compare our results with experimen- 
tal data on electron-doped cuprates. Note that at weak coupling the dominant correla- 
tions are strongly affected by the the Fermi surface topology and its intersection with 
the umklapp surface, while no qualitative changes occur at the transition from electron 
to hole-doped fillings. Hence, we may discuss qualitative properties of electron-doped 
cuprates, although we computed for densities n < 1. Unfortunately, the angular- resolved 
photemission (ARPES) data for electron-doped cuprates available to this day are con- 
troversial. Very recently Claesson et al.— have presented measurements which do not 
show any signs of pseudogap behaviour, or other fingerprints of strong correlations in the 
(tt, 7r)-transferring channel. On the other hand, earlier observations by Armitage et al.— 
are indeed partially consistent with our results. There, in a certain parameter range the 
Fermi surface vanishes around the hot spots where it would intersect the umklapp surface, 
and single-particle spectral weight is shifted to lower energies in that region. However, 
in these samples the qualitative picture changes as a function of doping. For a heavily 
underdoped system small electron pockets appear near (vr, 0), while the remaining part 
of the would-be Fermi surface is destroyed, also along the zone diagonal. Toyama and 
Maekawaii^ argue from calculations on a. t — t' — t" — J model that this may be due to the 
closeness to the antiferromagnetic phase, which is not as easily destroyed by doping as it 
is for hole-doped cuprates. Alternatively, this feature can be due to a small change in the 
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repulsive interaction as a function of doping, as suggested in reference 

the interaction needs to be sufficiently large, and is outside the weakly coupled regime. 

A comprehensive understanding of cuprate superconductors and other strongly corre- 
lated electron systems is still hampered by a lack of theoretical methods and frequently 
also by difficulties in obtaining reliable experimental data. At the present stage of its 
development the fRG can contribute by offering a controlled and transparent treatment 
of electron correlations at weak coupling. This captures already a variety of non-trivial 
phenomena, but misses genuine strong coupling physics, such as the Mott transition. 
In the long run one may envisage combining strong coupling methods for short range 
correlations with fRG schemes for the long range part. 
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FIG. 1: Flow equations for the self energy S and the two-particle vertex F , respectively; the 
internal lines without slash correspond to the bare propagator D , the lines with a slash to its 
A-derivative d\D . 




FIG. 2: Fermi surfaces (solid lines) for n = 0.92 and n = 1 in the first quarter of the Brillouin 
zone, Fermi points used for the discretization of momentum dependences, and the first 6 mo- 
mentum space patches (confined by broken lines from (0,0) to (7r,7r)). The umklapp surface 
(broken line from (0, tt) to (vr,0)) intersects the Fermi surface at the so-called hotspots. 
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FIG. 3: Imaginary part of the self energy IniS(ki?,tj) normalized by U'^ as a function of lo for 
n = 0.92 and k_p on the hot spot (patch 1, see Fig. 2). The result from the fRG (blue solid line) 
is compared to the result from second order perturbation theory (black dashed line). 
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FIG. 4: T,(kp,Lo) normalized by C/^ at low energies lu for n = 0.92 and k^? in patches 1, 2, . . . , 6. 
Two results from the fRG with F^j^^ = 150 (blue solid lines) and F^^^^ = 15 (red dashed lines) are 
compared to each other; the former is obtained by choosing U = 1.815, the latter for U = 1.685. 
The results from second order perturbation theory for U = 1.815 (black dash-dotted lines) are 
also shown for comparison. 
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FIG. 5: Spectral function A(kF,uj) at small frequencies iv as obtained from the self energy 
results presented in Fig. 4. 
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FIG. 6: T,{kF,u;) as in Fig. 4 but for n = 1. 
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FIG. 7: Spectral function as in Fig. 5 but for n = 1. 



